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Preface to First Yearly Report under Contract NAS2-7613 

Work under Contract NAS2-7613 started on July 1, 1973, 
It is a continuation of research conducted since February 
1, 1967, under Contract NAS2-4151, Phase VII-A, B and C 
Reports of June, 1973, titled, "Concepts for a Theoretical 
and Experimental Study of Lifting Rotor Random Loads and 
Vibrations", are the final reports under Contract NAS2-4151 
and list the 9 preceding reports and 11 published articles 
and papers prepared under the contract. 

Meanwhile 2 further papers generated under this Con- 
tract have been published: 

Crews, S. T., Hohenemser, K. H. and Ormiston, R. A., 

"An Unsteady Wake Model for a Hingeless Rotor", Journal of 
Aircraft Vol. 10, llo. 12, Dec. 1973, pp. 758-760. 

Hohenemser, K. H. and Prelewicz, D. A., "Computer Experi- 
ments on Periodic Systems Identification Using Rotor Blade 
Transient Flapping-Torsion Responses at High Advance Ratio", 
AHS/NASA Ames Specialists Meeting on Rotorcraft Dynamics, 
Moffett Field, California, February 1974. 

The first paper which has been generated under the new 
Contract NAS2-7613 is: 

Hohenemser, K. H. and Yin, S. K. , "On the Use of First 
Order Rotor Dynamics in Multiblade Coordinates", presented 
at the 30th Annual National Forum of the American Helicopter 
Society, May 1974, Preprint 831. 


II 

The research goals stated in Contract NAS2-7613 are 

(a) Assess analytically the effects of fuselage motions on 

V 

stability and random response. The problem is to 
develop an adequate but not overly complex flight 
dynamics analytical model and to study the effects 
of structural and electronic feedback, particularly 
for hingeless rotors. 

(b) Study by computer and hardware experiments the feasi- 
bility of adequate perturbation models from non-linear 
trim conditions. The problem is to extract an 
adequate linear perturbation model for the purpose of 
stability and random motion studies. The extraction is 
to be performed on the basis of transient responses 
obtained either by computed time histories or by 
model tests. 

(c) Extend the experimental methods to assess rotor wake- 
blade interactions by using a 4-bladed rotor model 
with the capability of progressing and regressing blade 
pitch excitation (cy''lic pitch stirring), by using a 
4-bladed rotor mode with hub tilt stirring, and by 
testing rotor models in sinusoidal up or side flow. 

The first yearly report under Contract NAS2-7613 is 

subdivided into 3 parts, whereby Parts I, II, and III are 

related to the research goals (a), (b), and (c) respectively. 

The authors and titles of the three parts are: 


i 

I 


I 


I 

I 
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Part I, Hohenemser, K. H. and Yin* S. K., "Methods Studies 
Toward Simplified Rotor-Body Dynamics". 

Part II, Hohenemser, K. H. and Yin, S. K. , "Computer 

Experiments in Preparation of System Identification 
from Transient Rotor Model Tests". 

Part III, Hohenemser, K. H. and Crews, S. T. , "Experiments 
with a Four-Bladed Cyclic Pitch Stirring Model 
Rotor" . 

Part I considers a number of simplifications in rotor-body 
dynamics and applies the various analytical models 
to a hypothetical compound hingeless rotorcraft 
with and without feedback into cyclic and collective 
controls . 

Part II deals with the problem of rotor parameter identi- 
fication from noise polluted transient blade flapping 
responses. Computer experiments are used in order 
to gain some insight into the efficiency of various 
identification schemes to be later applied to rotor 
model flapping transients. 

Part III summarizes the test results obtained with the 
4-bladed cyclic pitch stirring model rotor. The 
analytical blade flapping responses without con- 
sidering the rotor wake are compared to the measured 
responses which include the wake-blade interactions. 
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METHODS STUDIES TOWARD SIMPLIFIED ROTOR-BODY DYNAMICS 
Part I of First Yearly Report under Contract NAS2-7613 

Abstract 

This report is directed to the problem of developing an 
adequate but not overly complex linear :'light dynamics analytical 
model of a rotorcraft to study stability, control, gust and 
random turbulence responses. Since the conventional flight 
dynamics analysis using quasisteady rotor derivatives is ade- 
quate for the long period modes like the phugoid mode, only 
short time responses are considered here, where rotor-body 
coupling is of importance. Thus the body motion consists of 
pitch, roll and vertical motion, omitting linear longitudinal 
and lateral and yaw perturbations. Five analytical models of 
varying degree of sophistication are applied to a hypothetical 
hingeless compound helicopter operating up to .8 rotor advance 
ratio. Stability and response data are obtained for the basic 
helicopter and for the vehicle with two simple control feedback 
systems. The least sophisticated analytical model which produces 
adequate results is determined for the various tasks. The 
random vehicle response to atmospheric turbulence shows a 
remarkable attenuation with either of the two control feedback 


systems 
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Nomenclature 
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E[...] 

F 

I 

K 

L 


M 

N 

P 


R 

Z 


tip-loss factor 
expected value of ••• 
state matrix 
moment of inertia 
feedback constant 

right rolling moment over A/C rolling 
inertia, unit also turbulence 
scale length 

nose-up pitching moment over A/C 
pitching inertia, unit R^ 

blade number 

blade elastic first flap-bending 
frequency when rotating. unit R. also 
covariance matrix 

rotor radius 

down normal force over A/C mass, unit 
Rr2 


m 


n 

P 

q 

t 

w 


Y 


rotorcraft mass 
blade mass 

zero mean white noise 
right rate of roll, unit R 
nose-up rate of pitch, unit R 
time, unit 1/R 

down normal velocity, unit RR 

blade flapping angle defined by line 
from rotor center to blade tip 

blade Lock number referred to rotor 
center 




'.’'I 

,5 

k 





- ^ 
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6 

6 

X 

V 

C 

a 

T 

♦ 

* 

W 

n 

Superscripts ! 


T 


control input 

blade pitch angle 

up gust velocity, unit RQ 

advance ratio 

real part of root 

standard deviation 

time constant (time unit 1/Q) 

control phase angle 

blade azimuth angle 

« 

frequency, unit (I 

rotor angular speed, rad/sec 


time derivative, time unit 1/Q 
transposed matrix 


Subscripts : 

a augmented, or acceleration 

i integral feedback for rotor tilting 

X about longitudinal axis 

y about lateral axis 

I forward cyclic 

II left cyclic 


o 


collective 
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Part I 


Introduction 

In Reference [1] stability and short time responses of a 
winged hingeless rotorcraft are studied with a linearized 
system of equations which includes 3 rotor modes • coning, 
regressing and progressing - , which includes periodic 
coefficients and which includes first blade mode elastic 
deflections using the methods of References [2] and [3]. 

Three control feedback systems were considered: Coning feedback, 

proportional cyclic feedback and a combined feedback. It 
was found that all three feedback systems • particularly 
the last one - were effective in alleviating control over 
sensitivity, control and damping cross coupling, pitch 
divergence and gust sensitivity. 

The results are useful as trend studies but are only 
approximate because inplane blade dynamics, rotor wake dynamics 
and horozontal rotor forces have been omitted. The justification 
for omitting inplane blade dynamics is that inplane blade 
oscillations affect the body motions only in a narrow frequency 
band centered at the air resonance frequency. This fact 
evolved rather clearly in recently obtained unpublished 
frequency responses of the B0>105 helicopter, where the air 
resonance frequency is .350, which is at the upper end of 
the flight dynamically important frequency spectrum. Coupling 





of a high gain feedback system or of the human controller 
with the air resonance mode may under certain conditions 
destabilize this mode, or at least cause high in> < blade 
loads. Apart from this special problem the ef :ct of . .ade 
inplane dynamics on flying qualities of rotorcraft is expected 
to be usually small. 

The justification of omitting rotor wake dynamics is 
that an equivalent reduced Lock number can approximate these 
effects » at least in the low end of the frequency spectrum 
which is of main importance in flight dynamics, see for 
example Reference [4]. Finally, the justification for 
omitting horizontal rotor forces is that their moments with 
respect to the aircraft C.6. are much smaJler than the moments 
transmitted from the blades to the hub, at least for the 
rather high blade flapping frequencies assumed for the 
numerical examples. Such high flapping frequencies can be 
expected for unloaded slowed rotor operation in cruising of 
compound rotorcraft. For the short time responses the lon- 
gitudinal and lateral linear perturbation motions of the 
rotorcraft have been neglected anyway, so that horizontal 
forces are not needed in the analysis. 

The question studied in this report is whether simpli- 
fications in the rotorcraft representation of Reference [1] 
can be applied without introducing sizeable errcrs. The 
motivation to simplify the analytical model as much a.> feasible 
is not so much to save computer effort in establishing time 



histories or stability characteristics, but rather to obtain 
a better visibility for the essential pa smeter effects* 

Another motivation is related to the fact that system 
identification algorithms and random response analyses do 
tend to require substantially more computer effort for 
increasing order and for time variability of the system. 

While the basic compound rotorcraft treated in this report is 
the same as in Reference Cl], different feedback systems 
have been assumed in order to cover a wider range of such 
systems. This report is a more elaborate version of 
Reference [S], containing more analytical details and more 
numerical data. 

Basic Rotorcraft Description 

Though the same rotorcraft is used as an example • in 
Reference [1], a brief description is given here f< 
sake of completeness. The hingeless rotorcraft is .he 
compound type, cruising with reduced rotorspeed in an 
unloaded rotor condition. It is assumed that the rotorcraft 
performs a uniform forward motion and is restrained in yaw 
and side motion. The body is free to pitch, roll and move 
vertically. Thus we have added roll to the usual short time 
longitudinal flight dynamics because pitch and roll are eojpled 
through the rotor modes. The classical phugoid, dutch roll 
and spiral modes cannot occur with the assumed restraint. 
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Thes« modes are little modified by coupling with the rotor 
modes and reasonable approximations can be obtained with the 
conventional derivative analysis. 

Due to the reduced rotor speed the blade flapping frequency 
in the rotating reference system is assumed to be rather 
high* P s 1.2, This assi'tiption leads to a highly unstable 
basic rotorcraft and will impose severe demands on the 
selected control feedback systems. The stiff rotor also 
simplifies the analysis since in the first approximation inplane 
rotor forces need not be considered because of the relative 
smallness of their moments about the rotorcraft c.g. as 
compared to the hub moments. For advance ratio g s .8 flap* 
bending flexibility of the blades and reversed flow effects 
are taken into account following Reference [3]. A single 
elastic blade flap*bending mode is used. For advance ratio 
g s ,4 flap*bending flexibility and reversed flow effects are 
neglected. The effect of the rotor wake dynamics on rotor 
frequency responses in the linear range was studied in Reference 
[4] from which one can conclude that for the high frequencies 
of the coning and advancing tilting modes dynamic wake effects 
are small, while for the low frequency regressing tilting mode, 
the assumption of a quasi steady wake is a reasonable 
approximation. The effects of the latter can be approximated 
by a certain reduction in blade Lock number, see References [4], 
[10] and [11]. Values for the blade Lock number y in the 


t 
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numerical examples of this paper are 5 and 8, covering the 
range of current hingeless rotors. 

The wing area is 6% of rotor disk area. The wing lift 
slope is assumed to be 4.5» which would include rotor*wing 
interference effects. The horizontal tail area in excess of 
that required to neutralize the instability of the fuselage 
without rotor is 1.5% of rotor disK area. Tail moment arm is 
1.2R, tail lift slope including wing and rotor downwash effects 
is 1.8. Downwash lag effects on the tail are neglected. With 
these values and using the rotor radius as length unit and Q 
as angular velocity unit, one obtains for the body derivatives 


at .8 and .4 advance ratio 


U 

- .8 

u « . 

Wing roll damping 


s 

-.0200 

-.0100 

Tail pitch damping 

Mq 

s 

-.0075 

-.0038 

Tail pitch stability 

Mw 

s 

-.0063 

-.0032 

Wing plus tail normal damping 

2w 

5 

-.0288 

-.0144 

Tail normal force from pitch 
rate 

2q 

s 

-.0031 

-.0015 


Body and blade masses and inertia moments are related by: 

“ 5, I^/X^ • 75, mj^/m ~ .20, • .60 

Data not shown here have also been obtained for ly/l^ - 25 with 
essentially the same conclusions. The rotor is assumed tu be 
3>bladed. The control phase angle for the basic rotorcraft is 45^. 
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Feedback Systems 

Two feedback systems are considered. The first* has 
rotor tilting plus body angular rate feedback with a large 
actuator lag. It was described in Reference [6] and its 
performance is almost identical to that of the Lockheed gyro 
control AMCS described in Reference [7], though the mechanic 
zation is different. Wind tunnel model tests with a variant 
of this system using an isolated hingeless rotor are described 
in Reference [8], In the flight dynamics frequency range this 
feedback system emasculates the rotor with respect to tilting* 
so that angle of attach instability* control over sensitivity* 
control cross coupling* cross damping, and stick reversal in 
g maneuvers are avoided. The feedback equations in their 
linear form are 


+ 0|/t = -KiBj + 

COS 

♦ 

- «ii 

sin ^ 

t Kjq 

(1) 

+ 0jj/t = -K.Bjj ♦ 6j 

sin 


+ 6ii 

cos ^ 

. Kjjp 

(2) 


The values for the parameters used in the numerical examples 
are the following. The actuator time constant is t = 10, 
which is the time for 1.6 rotor revolutions. The rotor tilting 
gain is = .3, the control ph,ise angle is ♦ s 7.5°, the rate 


gyro gains are = Kjj = 1. 


I' 


The second system has normal body acceleration feedback 
into collective pitch. Without fuselage aerodynamics this 
is equivalent to coning feedback. According to Reference [9] 
this system proved to be effective in substantially ex- 
tending the speed range for acceptable handling, characteristics 
of a hingeless rotorcraft. The feedback equation is 

Oo ♦ T Go = K^(w - qu) (3) 

For the numerical examples the actuator time constant is 
T s .5, which is the time for .08 rotor revolutions. The 
gain is Kg = 10. In addition one needs the cyclic control 
relations 

Qj = cos (j> - sin ^ (4) 

. Gjj s 6j sin <1 + fijj cos ♦ (5) 

where the control phase angle is (|> = 25®. The control phase 
angles for the basic rotorcraft and for the craft with the 
two feedback systems were selected in such a way as to 
minimize control cross coupling at y * 5 over the entire 
flight range between 0 and .8 advance ratio. For the craft 
with tilting feedbacx. the tail size was reduced to 1% of 
rotor disk area, sin_e this feedbacK system is more effective 
than the normal acceleration feedback. Schematics of the 
two feedback systems are shown in Fig. 1 



10 


Rotor/Body Equations of Motion 

One form of th<i linearized rotor/body equations of motion 
are given in Reference Cl], Eqs. (23) to (35). Here a more 
developed form is given which is directly usable for programming. 
For advance ratio u ^ rigid straight blades are assumed, 
which are elastically hinged at the rotor center. The 
equations for this case can be given in exact form, if the 
small reverse flow effects are neglected. For advance ratio 
ii s .8 reversed flow and blade elasticity are not negligible. 

In this case a Fourier analysis of the periodic terms was 
performed retaining the first 7 terms of the Fourier series. 

The equations were transformed to non rotating body fixed 
coordinates in state variable form. The basic rotorcraft is 

a ninth order dynamic system, with the state variables p, q, w, 

• • • 

$ 1 , 8j, &II* With tilting feedback according to 

Eqs. (1) and (2) one obtains an eleventh order dynamic system 
including the state variables 6^, With normal acceleration 

feedback according to Eq. (3) one obtains a tenth order 
dynamic system including the state variable 0^. 

At low advance ratios (u < .5), where rigid straight blades 
are assumed and the small reverse flow effects can be neglected, 
the linearized rotor/body equations of motion, for a three-bladed 
rotor, in terms of multiblade coordinates are 


J 


11 

B ♦ (bS/ 0)6 + (b\w/ 12)B., ♦ (B^YW^/16 )Bt si» 3^ 

O O XX X 

-(B^YU^/16)8 jj cos 3t + Bjj - 1.5 (w - wq) - (B\y/12)p 
-(b\/6)w = (B**y/8 + B^yu^/8)0^ - [bSp/® - 
(B^YU^/16) sin 3t]9j - (B^YU^/16)6 jj cos 3t (6) 

$j ♦ (b“y/8)8j + (p^ - l)8j + (B^YW/6)6^ ♦ 28jj 

+ (B**y /8 ♦ B^yu^/ 16 )Bjj + [(B^yu/ 6)6 j - (B\y/12)6jj] 
cos 3t + C(B^YU^/8)6p + (b\u/ 12)8 j + (b\p/ 6)8 j j ]sin 3t 

-q - 2 p - [B^y /8 ♦ (B^YP/12)sin 3t]q 

+(B^YW/12)p cos 3t = -(B^Yy^/8)0^ cos 3t 

+ (B^YU/8)8 j cos 3t + [B**y /8 + B^yp ^/18 + (B^YU/8)sin 3t]6jj 

( 7 ) 

6ii + (b'*y/8)8jj + (p^ - DBjj + (B^yu/ 6)6^ - 2§j + (B^yp2/16 

-b‘*y/8)6j + [(b\u/6)Bj - (B\w/12)8jj]sin 3t - C(B^YM^/8)6^ 

+(B^Yu/ 12)8 j + (B^Yy/6)8jj]cos 3t - p + 2q + (B^YM/12)q cos 3t 

-[(B^y/ 8) - (B^Yu/12)sin 3t]p - (B^YM/*^)w = [B^YP /3 

-(B^YP ^/8 )sin 3t]0 - [B*^y /8 + 3B^YP^/16 - (B^YP/6)sin 3t]0. 

o ^ 

- (B^YP/6)0 jj cos 3t 


( 8 ) 


12 


F 


li 


p - LpP = 1.5ajj/I^)(l - P )Bii (9) 

q - M^q - H^w = 1.5 (I^/Iy)(l - P^)0j (10) 


W - wq - z w - 2 q - (m./m)[1.5 0 (B'^y/6)S 

^ w q ^ D o o 

♦ (B^YU/e)(Bjj - Sj - p) + (B^YW/ 8 ) 0 j ♦ (Byu^/ 8 ) 

(Bt sin 3t - Bjj cos 3t) - (B^y/‘*)w] » (mjj/m)[-(B^Y/6 
+ Bym^/4)0^ ♦ (B^ym/‘*)0j + (Byw^/8)(0jj cos 3t - 0j sin 3t)] (11) 


n 




0 

!1 

n 


At higher advance ratios* e.g. y = .8, where blade 
flexibility and reverse flow effects have been considered, the 
coefficients in the linearized rotor/body equations of motion 
depend on blade flapping mode shapes and flow regions, and can 
no longer be expressed explicitly in general forms in terms of 
advance ratio y and blade tip loss factor B. Therefore the 
equations used in the numerical examples for y = ,8 are not 
given here. 

System Simplifications 

The linear rotor perturbation equations (6) to (11) are of 

second order in S.,* 8 and include periodic terms. We are 

1 XX o 

considering 5 types of system modeling indicated in Table 1. 
Type 1 is the complete system including periodic terms, body 
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motions, multiblade accelerations and multiblade rates, type | 

2 is the rotor periodic system obtained by rigidly restraining ] 

the body. Type 3 is the rotor/body constant coefficient 
system of second order in each multiblade coordinate 
obtained by omitting the periodic terms. With full body 
freedom this has been also called the 9x9 system since 
there are 6 body and 3 multiblade degrees of freedom. Type 4 
is new and is introduced in this paper and in Reference [5]. 1 

The system is of first order in each multiblade coordinate. I 

The periodic terms and the multiblade accelerations are omitted. i 

Type 5 is of zero order in the multiblade coordinates; both | 

multiblade accelerations and rates are omitted. This system 
is equivalent to the conventional 6x6 quasisteady derivative 
system. ] 

i 

1 

Table 1 | 

rive Types of System Modeling \ 



Periodic 

Body 

Accelerations 

Rates 

I 

Type No. 

Terms 

Motions 

•« 

6l 

*i «:i 8. 

1 


F 
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The effects of omitting the multiblade accelerations will 
be shown for the case of an isolated rotor in vacuum. The 
equations for the multiblade tilting coordinates Sj, Bjj are 

Bj + (P2 - 1)B^ ♦ 2 Bii = 0 (12) 

♦ (P2 - DBii - 2 Bi = 0 (13) 

which are satisfied for 

, Bjl = i * u = 1 1 P. (1*») 

The solutions describe an advancing tilting mode with 
frequency 1 + P and a regressing tilting mode with frequency 
P - 1. 

•« •* 

Omitting the accelerations Bj , Bji , the equations have 
only one root, u = (p2 - l)/2 which is approximately equal to 

P - 1 for not too large P - 1. Thus first order rotor 

dynamics in multiblade coordinates removes the advancing 
tilting mode with the high frequency 1 + P and yields a 
regressing tilting mode with approximately the same frequency 
as for the complete system of equations. 

For the coning node in hovering we have 

6,5 (y/BHo ♦ P^Bc. = 0 (15) 

Leavim^ out the r^ccelera t i on 6 reduces the coning degree 

of freedom to a first order system with spring and damping 

only. Note that what is called here first order rotor dy- 

«• 

namics is no t equivalent to omitting the S term in the 
individual biade flapping equation. The terms omitted are 
the multiblade accelerations 6,, Bjj an.. Bq . Physically 



We will consider here both the rotor/body periodic system 


and the isolated rotor periodic system. In Reference [1] it 
was shown that rotor/body coupling had little effect on the 
root curves for coning and proportional tilting feedback. 

Here the comparison is extended to the lagged tilting moment 
feedback expressed in Equations (1) and (2). The analysis 
makes use of the Floquet state transition matrix in the 
multiblade coordinate form and extracts the characteristic 
values from this matrix. The ambivalence in assigning frequency 
values makes it possible to use only the positive frequency 
half plane. For the 3-bladed rotor assumed in this study the 
flapping instability occurs as a coupled advancing and coning 
mode with frequency 1.5. Figure 2 shows the root plots for 
the rotor tilting feedback system with gain as parameter. 

The advance ratio is .8 and the cases of blade Lock number of 
5 and 0 are shown. There is very little difference in the 
root curves for the rotor alone - dashed lines - and the 
rotor/body system, solid lines. However, the blade Lock number 
has a substantial effect* whereby the limiting gain is reduced 
from 5 to 8. From Figure 2 it appears that a stability 
analysis for the coupled advancing and coning mode can be 

I 
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mad* with th* rotor alon* and rotor/body coupling will not 
affact th* r*sult very much. 

For normal acceleration feedback the case of the rotor 
alone is not really comparable to that for the complete 
rotor/body system, since the fixed wing contributes sub- 
stantially to the vertical damping thus improving the 
stability. Therefore, Fig. 3 shows only the case with 
rotor/body coupling for an advance ratio of .8. The 
instability of the coupled advancing and coning mode occurs 
in the same way as before with a frequency of 1.5. Increasing 
the blade Lock number from 5 to 8 reduces the limiting gain 
from about K = **0 to 30. The stability of the coupled 
advancing and coning mode increases with decreasing advance 
ratio, where the coefficients of the periodic terms are 
smaller. The highest advance ratio is, therefore, the most 
critical one. 

Eigenvalue Comparisons 

Table 2 shows in the first column the eigenvalues for 
advance ratio u = .6 and blade Lock number y ® 5. It shows 
in the second coiumn the effects on eigenvalues of changing 
at Y = 5 the advance ratio from .8 to .4, and in the third 
column the effects ot changing at u = .6 the blade Lock 
number from y - i> tu 8. The complete equations of the 
periodic system were used. The fourth column gives for 
U » .8 anu Y = S the effects of omitting tne periodic terms 
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In the multiblade equations^ the fifth column ehowa the effects 
of omitting in addition the accelerations, and the 6th column shows 
the effects of omitting in addition the rates of the multi- 
blade coordinates. The 3 sets of rows in Table 2 refer to 
the 3 cases of the basic rotorcraft, of the craft with rotor 
tilting, and of the craft with normal acceleration feedback. 

The coupled modes are named according to the predominant 
uncoupled modal content of each. As mentioned before, the 
basic rotor/body system is of 9th order, with tilting feed- 
back of 11th order, with normal accsleration feedback of 10th 
order. With first order rotor dynamics in each multiblade 
coordinate the order of all systems is reduced by 3. With 
zero order rotor dynamics, representing the derivative approach, 
the order is reduced by 6, 

To look into the effect of advance ratio, compare the 
first two columns of Table 2, The basic rotorcraft has at 
u s . 8 a pitch (iiverf’ence which is about 6 times more 
severe than at u = . *• . The other eigenvalues change relatively 
little with ii • v;ith tilting feedback the pitch divergence 
is replaced by a stable oscillation. One roll mode is 
more damped, the other less clamped than without feedback. 

Coning and advancing modes show little change with feedback 
or with advance ratio. With normal acceleration feedback 
the pitch divergence is also removed, however here the pitch 
mode changes from two convergences to an oscillatory mode 
when going from u = .'♦ to u = .8, 
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To look into the effect of advance ratio* conpare the 
firat two columns of Table 2« The basic rotoreraft has at 
w ■ .8 a pitch divergence which la about 6 times more 
severe than at w ■ *4. The other eigenvalues change relatively 
little with u. With tilting feedback the pitch divergence is 
replaced by a stable oscillation. One voll mode is more 
damped* the other less damped than without feedback. Coning 
and advancing modes show little chan, >m 'ith feedback or with 
advance ratio. With normal accelerat * feedback the pitch 
divergence is also removed* however here the pitch mode 
changes from two convergences to an oscillatory mode when 
going from u « .4 to u « .6. 

^ To look into the effect of blade Lock number at v ■ *8 
compare the first and third columns of Table 2. The basic 
rotoreraft shows a large increase in pitch divergence from 
increased y. With tilting feedback the blade Lock number 
has little effect on the pitch mode* while some of the other 
modes are strongly affected by y. The same is true for the 
craft with normal acceleration feedback. 

To look into the degree of approximation provided by 
the 3 constant system models* compsre the first column with 
the last 3 columns of Table 2. The pitch mode eigenvalue is 
in most cases reasenabl/ well approximated even by the sero 
order system. The two roll mode eigenvalues are well 
approximated by the first order system in each multiblada 


coordinate, but considerably in error for the zero order 
system. The second order syscem provides reasonably 
accurate eigenvalues for all modes, however data not shown 
here indicate that with increasing blade Lock number 
Increasing effects of periodicity occur. 

Step Control Response Comparisons 

All transient responses shown in Figs. 4 to 12 are given 
for a non-dimensional time period of 40 which is the time 
for 6.3 rotor revolutions. One can assu;..e that within this 
short time period the effects of the omitted phugoid, dutch 
roll and spiral modes (7ill be relatively small as compared 
to the effects of the short period pitch and roll modes. 

The figures show 4 response variables: The roll rate p, the 

pitch rate q, the normal velocity w, which in a body fixed 
reference system is proportional to angle of attac , and the 
normal acceleration w - qu. In Figs. 4 to 10 the responses to 
both unit longitudinal and lateral control input are shown. 

Fig. 4 is for the basic rotorcraft at y » .8, y * 5, 

# * 45®. The high degree of instability and the high control 
cross coupling is evident in spite of an effort to minimize 
the control cross coupling by selecting a proper control 
phase angle y. The oscillations in the roll rate and 
normal acceleration responses have a frequency of 3 and are 
caused by the periodic terms. In a linear analysis these 
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3 per rev. oscillations for thd 3-bladed rotor must be 
superimposed to the oscillations from the trim condition. 

Fig. S, also for the basic rotorcraft at the same conditions, 
shows the responses to the unit control inputs for the 
constant first order and zero order systems. The first 
order system responses are very close to those of the 
periodic system in Fig. 4, except for the missing oscil- 
lations. Use of the zero order system leads to substantial 
errors. 

Fig. 6 shows for the rotorcraft with tilting feedback 
at y s .8, Y * 5, ♦ = 7.5® the responses to unit control 
inputs computed with the zero and first order systems. 

Pitch and roll rates stabilize within a short time period, 
the normal acceleration soon becomes concave downward. 

Control cross couplings are negligible. The zero order 
responses deviate initially from the first order responses, 
particularly in roll rate. For the remaining numerical 
examples the first order approximation was used. Fig. 7 shows 
for the rotorcraft with tilting feedback the effects of blade 
Lock number on the responses to unit control inputs for 
y s ,8, t * 7.5°. It is seen that this feedback system is 
not sensitive to changes in blade Lock number which agrees 
with the concept of "rotor emasculation". 


J. 
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rig. 8 shows for the rotorcraft with normal acceleration 
feedback the effect of blade Lock number on the responses to 
unit control inputs for u - *6, ^ « 25^* The improvement 
ovar the basic rotorcraft is very good for y ^ S. Pitch and 
roll rate stabilize rapidly* and the normal acceleration 
becomes concave downward after about 4 rotor revolutions. 
Control cross coupling is small. For y - 8* the roll rate 
from longitudinal control input is large* indicating that 
^ s 25® is not optimal. 

The effects of advance ratio on the step control responses 
are shown in Figs. 9 and 10. Advance ratios v * *4 and .8 
are assumed. Fig. 9 is for the rotor tilting moment feedback. 
It is seen that control cross coupling remains small for 
. both advance ratios. The normal acceleration per unit control 
input is for u = .4 smaller than for y * .8. Fig, 10 is for 
the nox'mal acceleration feedback. Here the case for y « .4 
shows more control cross coupling than the case for y - .8* 
indicating that for this type of feedback the selection of a 
control phase angle which keeps the control cross coupling 
effects small over the entire flight regime* is somewhat of 
a problem which does not exist for the ”emasculated” rotor 
with the rotor tilting feedback. 

Step and Random Cust Response Comparisons 

Gust responses are obtained from the equations of motions 
by replacing the normal velocity w by the sum of normal 
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body velocity and normal gust velocity w f X, whereby X 
is prescribed. For unit gust step inputs X is equal to 
the unit step function. For random gust input the standard 
deviation O)^ is assumed to be a unit step function of time. 
The standard deviations for roll* pitch, normal body 
velocity and normal body acceleration, Op, Oq, Ow, 
respectively are then computed for the case of the rotorcraft 
encountering at time t = 0 a turbulence region with unit 
standard deviation Since the basic rotorcraft is highly 

unstable, a steady turbulence response state does not 
exist, however, the initial random responses can be compared 
to those with control feedback. The gust responses shown 

f 

in Figs. 11 and 12 have again been obtained with first 
order rotor dynamics 

The random response analysis follows Reference Cl2j. 

The "point" assumption is used, according to which the 
entire rotor disk experiences at a given time the turbulence 
velocity at the rotor center. The normal gust velocity is 
obtained by passing white noise through a first order filter 

X + aX = ox(2a)^^^n (16) 


where 


a s 2u/(L/k) 


(17) 
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Th« ratio of turbulence scale length L over rotor radius R 
is assumed to be 12, for example L s 350 ft, R » 3C f\. The 
rotor/body equations of motion written in multiblade state 
variable form are augmented by Cq. (16) with the additional 
state variable X. The augmented state vector Xq has the 
components 0^, 0^, P* random gust response 

equations then read 

Xg - FaXg ♦ G V (18) 

where Fg is the augmented state matrix and where G and v are 
given by 


0 • 

• • 0 


m * 

n 

• 

• 


• 

• 

« 

V » 

• 

• 

0 • 

. . i2a)i/2 


• 

n 

m • 


n is zero mean white noise with standard deviation 0 |j * 1* 

The equation for the augmented covariance matrix Pg is, see 
for example Reference [12], 

P =FP tPF’’ + GG^ ^20) 

d a a a a 
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For purposes of detern Inin? the response standard deviations 
we need only the unaurmented 6x6 covariance matrix without X 

P * Rxx * KLxx^J (21) 

and the unaurnent *‘d .state vector equation 

X = Fk (22) 

The standard deviatir>ns for th*.* stfte vario' o3 Op , Oq, o^, 
are directly obtained from the diagonal terms of P * Rxx* 

The covariance matrix for the normal acceleration is 

Rgg s Et (w-qy ) (w-qu )'^3 

= Kw - ‘W^wq ♦ W^Rqq (23) 

Rq*q is a component of the covariance matrix Rxx« Rwq is s 
component of 

RAx * r Rxx (2*» ) 

and is a component of 

Rxx = RxxF^ (25 ) 

Thus the standard deviation of the normal acceleration 
can be computed from Eqa. (23) to (25). 

Fig. 11 shows the unit step gust responses p, q, w, 
w - qw for y a .0 and y * 5. Responses for the basic rotor- 
craft and for the craft with the control feedback systems are 
given. In addition, the q, w, and w - cm responses are also 
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shown without the rotor. The combination of body, wing 
and horizontal tail represents a stable configuration. 

The destabilizing effect of the rotor is clearly recog- 
nizable, also the stabilizing effects of the two feedback 
systems. The initial maximum normal acceleration is not 
much affected by the feedback systems, since a large 
portion of this normal acceleration has its origin in the 
fixed wing. 

rig. 12 shows the standard deviations of the same 4 
variables following entry into a turbulence region at t * 0. 
The feedback systemr are very effective in reducing the 
random responses in roll and pitch rate. Without rotor 
the pitch rate standard deviation Oq is almost as high as 
with rotor. Either of the two feedback systems brings a 
large reduction in random pitch rate response even as 
compared to the stable configuration without rotor. The 
normal acceleration standard deviation Ow-qy is also caused 
in part by tne fixed wing. The feedback systems reduce the 
rotor contribution considerably, but the overall value is 
not decreased by a large amount. In alleviating normal 
accelerations from cither step or random gusts, the normal 
acceleration feedback system is slightlv more effective 
than the rotor tilting feedback system. The normal velocity 
standard deviation is highest without rotor and lowest for 
the craft without ''cedback. T>iis is a consequence of the 
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rotor normal dampin?. With feedback systems rotor r.ormal 
dampinf» is reduced and intermediate values of normal 
velocity standard deviation obtained. 

Conclus i ons 

The following conclusions are based on results of a 
linear coupled rotor/body flight dynamics analysis for a 
hingeless compound helicopter operating up to .8 advance 
ratio. The hingeless rotor is rather stiff in flapping 
(P s 1.2), the fixed wing has 6% rotor disk area, the 
horizontal tail has 1.5% rotor disk area in excess of that 
required to neutralize the body instability. The conclusions 
are subdivided into those on analytical modeling and those on 
flying characteristics of the hypothetical rotorcraf t . 

Conclus i ons on Analytical Modeling 

• The analysis of flapping instabilities due to inter blade 
coupling from the control feedback systems requires the 
retention of the terms with periodic coefficients in the 
equations of motion. For rotor tilting feedback into 
cyclic pitch, rotor /body coupling effects are small and 
the isolated rotor analysis is adequate. For normal 
acceleration feedback into collective pitch, damping of the 
vertical motion by the fixed wing is beneficial and a 
coupled rotor/hody flapping stability analysis is required. 


26 


•• The analysis of transient responses to cyclic pitch 

step inputs including the terms with periodic coefficients* 
resulted In substantial N per rev. oscillations of the 
rigid body In roll rate and in normal acceleration. If 
these rigid body oscillations are of no Interest* constant 
coefficient modeling can be used. Of the three constant 
coefficient models studied* the one using second order 
rotor dynamics In each multiblade coordinate • equivalent 
to what is sometimes called the 9x9 model - was found to 
be over-sophisticated, the one using first order rotor 
dynamics was found to be adequate, the one using zero order 
rotor dynamics - equivalent to the 6 x 6 or derivative 
mov^el - can result in substantial errors. If the rotorcraft 
is stabilized by a- feedback system, these errors are 
limited to the initial response and fade out with time. 

For an unstable craft the errors increase with time. 

Conclusions on Flyinp Characteristics 

• Following a cyclic pitch control step input the basic 

rotorcraft shows rapid divergence in roll rate, pitch rate* 
normal velocity, and normal acceleration. Large control 
cross coupling effects are unavoidable inspite of optimization 
of the control phase angle. Of the eigenvalues only that 
for a predominant pitch mode indicates divergence. Due to 
cross coupling a lateral control step input also leads to 
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rapid d£vargence» A vertical step gust leads to pitch 
divergence associated with considerable roll* Entering a 
region of vertical turbulence the standard deviations for 
roll rate, pitch rate, normal velocity and normal accel- 
eration rapidly reach a temporaty level from which they 
then gradually increase. About 60% of the normal accel- 
eration turbulence response is caused by the fixed wing* 

• The rotor tilting feedback into cyclic pitch, even with a 
reduced tail area of 1%, completely stabilizes the rotor- 
craft and removes the cross control coupling effects for 
both Y s 5 and 8 and for u s .4 and .8* In a step gust 
this feedback does not alleviate the first normal accel- 

. eration peak, however the normal acceleration soon approaches 
zero rather than increasing again due to pitch-up as for 
the basic craft* When penetrating a turbulence area the 
final level of standard deviations in roll and pitch rate 
are only 30 to 40% of the initial level for the basic craft* 
The standard deviation level for the normal acceleration is 
not reduced materially since the fixed wing provides the 
major contribution* 

• The normal acceleration into collective pitch using now 
1*5% tail area also completely stabilizes the rotorcraft* 

For blade Lock number y S and advance ratio u s .6 cross 
control coupling effects are negligible, for y * B and 
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)i « .S, or for Y * 5 and n « ,4 using tha tarns control 
phase angle ther is substantial roll froa longitudinal 
control input. Otherwise this feedback system provides 
flying characteristics very similar to those for the rotor 
tilting feedback system and leads to somewhat lower values 
of turbulence response in roll and in normal acceleration. 
It is rather remarkable that for both feedback systems the 
pitch rate turbulence response is lower than without rotor. 
The lifting rotor with its feedback system absorbs part of 
the pitch rate turbulence response the craft would 
experience without rotor. 
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Figure Captions 

Sehtmatics of the Feedback Systems 

Root Plots for Rotor Tilting Feedback, ii ■ .8, Y ■ 5 
and 8, Periodic System 

Root Plots for Normal Acceleration Feedback, w ■ .8, 

Y « 5 and 8, Periodic System 

Control Responses of Basic Rotorcraft, |i « .8, Y * St 
8 ■ 45®, Periodic System 

Control Responses of Basic Rotorcraft, u ■ .8, Y * 5, 

8 « 45®, First and Zero Oraer Systems 

Control Responses with Rotor Tilting Feedback, w > .8, 

Y a 5, ♦ s 7.5®, First and Zero Order Systems 
Effect of Blade Lock Number on Control Responses with 
Rotor Tilting Feedback, 8 * 7.5®, u ■ .8 

Effect of Blade Lock Numoer on Control Responses with 
Normal Acceleration Feedback, 8 ■ 25®, w * .8 
Effect of Advance Ratio on Control Responses with 
Rotor Tilting Feedback, y • 5, 8 * 7.5® 

Effect of Advance Ratio on Control Responses with 
Normal Acceleration Feedback, Y * 5, 8 ■ 25® 

Step Gust Responses, v * .8, Y * 5, First Order System 
Turbulence Responses, u * .8, y * 5, First Order System 
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